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ABSTRACT 

A resonator  formed  by  nonperiodic  Bragg  structure  and  nonlinear  film  is  considered. 
Effects  of  the  second  or  the  third  harmonic  generation  are  theoretically  investigated. 
System  of  nonlinear  equations  for  amplitudes  of  time  harmonics  is  obtained. 
Phenomenon  of  resonant  increasing  of  amplitude  of  the  higher  harmonic  caused  by 
combination  of  nonlinear  and  Bragg  effects  is  discovered.  Reason  of  instability  of 
numeric  algorithm  for  “optimal”  case  is  revealed.  Ways  for  overcoming  this  instability 
are  found.  Bistable  behavior  of  scattered  waves  is  observed  for  cubic  nonlinearity. 

INTRODUCTION 

If  a harmonic  plane  wave  falls  on  a nonlinear  film  scattered  field  contains  higher 
harmonics  (the  second  harmonic  for  quadratic  non-linearity  or  the  third  one  for  cubic 
non-linearity).  As  a rule  values  of  non-linearity  and  amplitude  of  fundamental  harmonic 
is  not  very  large  so  that  amplitudes  of  generating  higher  harmonics  are  usually  small 
[1].  Higher-harmonic  amplitude  is  increased  under  phase-matching  condition  that  can 
be  achieved  more  easily  in  nonlinear  spatially  periodic  waveguiding  structure  [2], 
Amplitude  of  wave  falling  onto  nonlinear  film  can  be  increased  also  by  situating 
nonlinear  film  into  resonator.  Non-periodic  Bragg  structures  used  as  resonator  reflector 
[3]  have  flexible  frequency  characteristics  so  that  they  can  be  useful  for  devices 
operating  on  some  frequencies  simultaneously. 

METHOD  OF  ANALYSIS 

In  the  paper  physical  effects  in  a resonator  formed  by  non-periodic  Bragg  reflector  and 
non-linear  film  situated  on  metal  plane  are  theoretically  considered.  The  boundary-value 
problem  for  Hx  contains  spatial  inhomogeneity  in  wave  equation  and  non-linearity  in 
quasiimpedance  boundary  condition 


d2H  IdzdH  e(y)d2H 
dy 2 £ dy  8y  c 2 dt 2 
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where  wo  is  unperturbed  values  of  impedance  parameter,  wN  is  parameter  describing 
non-linearity  of  film,  q = 2 (quadratic  non-linearity)  or  q = 3 (cubic  non-linearity).  For 
quasiperiodic  Bragg  structure  permittivity  as  function  of  coordinate  y has  form: 

s(y)  = Efl  + p^£./(pT)exp[/Vy(y)J  (0 <y<yb)  (3) 

j 

T/t)  = Ix/(Pt)4f>  y>>  is  coordinate  of  Bragg-structure  beginning,  6 is  small  parameter  [3]. 
If  the  non-linearity  is  very  small,  undepleted  pump  approximation  could  be  used  [4]. 
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This  approach  is  easy  for  implementation,  but  it  does  not  allow  one  to  take  into  account 
changing  energy  of  fundamental  harmonic.  We  suppose  that  plane  time-periodic  wave 
(T  = 27r/co  is  time  period)  falls  on  the  structure.  Along  with  real  function  H(y,t)  we 
consider  complex  function  Hc(y>,t).  The  solution  of  the  boundary-value  problem  (1).  (2) 
for  H(y,t ) and  Hc(y,t ) wall  be  searched  as  a sum  of  monochromatic  w'aves 

N 

Hc{yj)=  X"cv,(T)exp(mcot).  m„(y)  = arn(y)  + iam(y\  (4) 

/;=! 


i\ 

H(y,t)  = Re[//,(y,/)]  = Xk,(T)cos(«ro/)-«,„(v)sin(wo/)], 


(5) 


n=\ 


A linear  spatial  problem  must  be  solved  for  each  frequency  harmonic  to  obtain 
transmission  matrices  of  Bragg  structure. 
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If  y>yb  functions  acn(y)  express  propagation  of  direct  and  opposite  plane  w'aves 

aJy): 


(6) 


«/,,/„  ex 


pKt  ~ T,  )kh  ] + ahnn  exp[-  i(  v - y,  )kh  J 


where  kh  = k0  e (yh ) , um,>  and  are  complex  amplitudes  of  direct  and  opposite  plane 


waves  in  beginning  of  Bragg-structure  (y  = y/,). 

Result  of  solving  equations  (6)  is  transmission  matrices  T„  , that  couple  amplitudes  of 
direct  and  opposite  waves  in  beginning  and  end  of  Bragg  structure 


where  ae(j„  and  ae0„  are  amplitudes  of  direct  and  opposite  waves  on  surface  of  nonlinear 


film  (y  = 0).  They  are  connected  with  values  acn o = ac„( 0)  and 

dy 

following  relations: 


T = 0 


by 


flcnO  gl 'tin  p clean . d cl,o  ikc(aet/n  acon ),  kc  ~ k^-jsjo)  . (8) 


Granting  nonlinearity  of  boundary  condition  it  is  more  convenient  to  operate  with  real 
function  H(y,t)  and  real  spatial  functions  an,(y)  and  a ,„(}’)  in  (2).  Taking  into  account 
expression  (5).  boundary  condition  (2)  can  be  written  in  the  form 

F(/)S  Hd -k2w0H„ -k2w„qHH  % -l)H2  + H0H„\=Q,  (9) 

where  H&  //o,  Hh  Hu  are  auxiliary  functions  of  time  coordinate  / 

N .V 

H,i  = Eh'o  n cos(«©/)-  a'l0n  sin  (mo/)].  H0  = £ko„  cos(wo/)-  ai0ll  sin(wo/)J 

H = l 11  = 1 
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H,  = “S"ko»  sin (mo/)+  a, 0„  Costco/)}  H„  = - jV  ko»  cos(wo/)-  al0n  sin(»w/)} 
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Real  parameters  a V„o.  ar„ o,  a ’,„o.  am o are  expressed  in  terms  of  complex  parameters  a ’cM 
and  acn 0 

ri/v/O  Rc( a(  „o).  <7/„o  I m(k /,7io).  c/  Rc(r/  a? o).  r/  m\)  1 m(  a end)* 

Introducing  frequency-harmonic  expansion  into  boundary  condition  (9)  in  discrete  time 
points  on  time  period  we  get  system  of  non-linear  equations  w'ith  respect  to  harmonic 


AV/ii*.  Ukraine,  IX-tii  International  Coneerence  on  Mathem-iticai.  Methods  in  Ei.e<  tromaoneth  ■ Theory 


MMET*02  Proceedings 


553 


amplitudes  on  output  from  Bragg  structure.  Here  we  consider  that  amplitude  of  falling 
fundamental  wave  is  specified  on  outer  interface  of  Bragg  reflector.  Obtained  non-linear 
system  has  been  solved  by  numerical  method  based  on  Newton  method  with  finite- 
deference  approximation  of  derivatives. 

THE  PHYSICAL  EFFECTS 

The  most  interesting  phenomena  occur  at  resonance  and  middle  values  (close  to  optimal 
ones)  of  non-linearity.  In  this  case  amplitude  of  fundamental  harmonic  near  the  film 
multiply  exceeds  the  amplitude  of  falling  wave,  the  amplitude  of  the  reflected 
fundamental  wave  on  outer  interface  of  the  resonator  is  close  to  zero,  energy  of  falling 
wave  almost  completely  transforms  to  energy  of  higher  harmonics.  Unfortunately  near 
optimal  values  of  parameters  quasi-Newton  method  collapses.  To  overcome  difficulty 
of  algorithm  convergence  we  describe  the  fundamental  harmonic  not  by  the  wave 
amplitudes  on  outer  boundary  of  the  Bragg  structure  but  by  the  amplitudes  of  the  direct 
and  reflected  waves  on  the  surface  of  the  non-linear  film. 

Obtained  solutions  essentially  depend  on  from  type  of  non-linearity.  For  cubic  non- 
linearity the  amplitude  of  the  wave  falling  onto  the  resonator  as  function  of  the 
amplitude  of  fundamental  harmonic  on  the  film  is  not  monotonic  thus  dependences  of 
amplitudes  of  all  scattered  waves  have  hysteretic  character. 

CONCLUSIONS 

Physical  effects  of  scattering  of  a monochromatic  plane  wave  on  non-linear  films 
situated  in  Bragg  resonator  are  theoretically  considered.  Mathematical  model  of 
processes  of  higher-harmonic  generation  are  derived.  Influence  of  phase  and  amplitude 
parameters  of  Bragg  reflection  on  amplitudes  of  scattered  waves  is  investigated. 
Conditions  of  appearance  of  bistable  regime  are  found.  Obtained  results  can  be  used  for 
designing  generator,  transformer  and  digital  devices. 
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